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Abstract. Radiative nucleon-deuteron capture and photo disintegration of the 
three-nucleon bound state with two-body final states are described. The de- 
scription uses nucleon degrees of freedom extended to include the excitation 
of a single nucleon to a A-isobar. The baryonic interaction and the electro- 
magnetic current couple nucleonic states and states with a A-isobar. Exact 
solutions of three-particle scattering equations are employed for the initial or 
final states of the reactions. The current has one-baryon and two-baryon con- 
tributions. The role of the A-isobar in the description of the considered photo 
reactions is discussed and found to be moderate. The spin observables A yy and 
T20 at 90° lab scattering angle can be calculated model-independently from the 
El Siegert term in the long-wavelength limit. 
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1 Introduction 

Trinuclcon photo reactions are theoretically described. Only reactions with two- 
body initial or final states are discussed, not three-body photo disintegration of 
the three-nucleon bound state. The considered reaction energies are well above 
three-nucleon break-up, but remain otherwise low, i.e., below 50 MeV in the 
three-nucleon cm. system. 

The present calculation of trinucleon photo reactions extends the descrip- 
tion of elastic and inelastic nuclcon-dcutcron scattering given by some of the 
authors in refs. [1-4]. Nucleon-deuteron scattering is studied in general [5], 
in order to search for unambiguous signals of a three-nucleon force. Indeed, 
discrepancies between theoretical predictions, derived from two-nucleon poten- 
tials, and experimental data got uncovered [3, 6]. In the same spirit trinucleon 
photo reactions are studied, in order to find clear indications for the working 
of a three-nucleon force and of corresponding two- and three-nucleon electro- 
magnetic (e.m.) exchange currents. 

The theoretical description of trinucleon photo reactions, given in this pa- 
per, is realistic in the following sense: A realistic nuclear interaction is employed 
when solving for the hadronic states involved, i.e., for the trinucleon bound 
state and for the nucleon-deuteron scattering states. The employed nuclear in- 
teraction is a hermitian coupled-channel two-baryon potential between nucleons 
which allows for the excitation of a single nucleon to a A-isobar. The A-isobar 
mediates an effective three-nucleon force. However, the Coulomb interaction 
between the two protons is not taken into account. Thus, the calculations of 
this paper are done for photo reactions on 3 H, though most experimental data 
will refer to photo reactions on 3 He. The nuclear current is built-up from one- 
baryon and two-baryon contributions. The current contributions contain, for 
reasons of consistency, pieces with transitions to one-baryon and two-baryon 
states with a single A-isobar. The description of the considered reactions is 
comparatively transparent: As long as only a few multipolcs of the current de- 
termine the observables, only few total angular momenta are picked out from 
the nucleon-deuteron scattering states. 

The calculation of trinucleon elastic and inelastic nucleon-deuteron scatter- 
ing, given in refs. [1-4], is based on a separable expansion of the two-baryon 
transition matrix; the same separable expansion is used here for the descrip- 
tion of photo reactions. In the separable treatment of the hadronic interaction, 
our calculation corresponds to the work of refs. [7, 8]; however, it is more com- 
plete with respect to the employed current. Our calculation attempts to be as 
realistic as the description given in refs. [9, 10], which does not require a sep- 
arable expansion of the hadronic interaction. In contrast, the addition of the 
A-isobar degree of freedom gives our calculation a new theoretical dimension. 
The calculation will attempt to isolate the effects arising from the A-isobar. 
Preliminary results were presented in ref. [11]. 

Sect. 2 describes the calculational apparatus. Sect. 3 shows and discusses 
results. Sect. 4 is a summary with conclusions. 
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2 Calculation Procedure 

2.1 Basic Elements of Calculation 

The basic building block for the description of photo reactions is the e.m. 
hamiltonian 

Ht m - = ljd 3 x J^(x)A fl (x)\ Xo=0 , (2.1) 

which couples baryonic states to the photon (7). The operators in Eq. (2.1) 
depend on space-time x, but are to be used as Schrodinger operators at time 
x = ; A^x) is the photon field operator, which we parametrize in the form 

(inY^hc f rf 3 fc 7 A 

M*) -7^)372-7 T^t^o 

[a A (k 7 Mk 7 A) e -* fe ^ + 4(k 7 )e;(k 7 A)e^] 1^^, . (2.2) 

The polarizations A = 1,2 correspond to those of a real transverse photon; 
e M (k 7 A) are the polarization vectors, constrained by fc^e^(k 7 A) = 0. A single- 
photon state of definite momentum k 7 and polarization A, ^-function normal- 
ized, is |k 7 A) = a^(k 7 )|0). 

The e.m. current operator J^{x) acts in the baryonic Hilbcrt space which 
we take to have two sectors, i.e., a purely nucleonic one and one in which one 
nucleon (N) is turned into a A-isobar. The two Hilbert sectors are displayed 
in Fig. 1 for a three-baryon system as 3 He. We use the current operator in its 
Fourier-transformed form, i.e., 

J"(Q) = J d 3 x e* Q - x J"(aO| Xo = 0) (2.3) 

and employ — specializing to a three-baryon system — a momentum-space 
representation, based on the Jacobi momenta ( pqK ) of three particles in the 
definition of rcf. [1], i.e., 

(K'p'qV M (Q)|Kpq) =(5(K'-Q-K)(p'q'|.? Ai (Q,K' + K)|pq>. (2.4) 

In Eq. (2.4) Q is the three-momentum transfer by the photon; it will take on 
particular values depending on the considered reaction; in e.m. reactions with 
real photons it is given by the photon momentum k 7 . A total- momentum con- 
serving (5-function is split off; the remaining current operator j' i (Q,K / + K) 
only acts on the internal momenta of the three-baryon system with a paramet- 
ric dependence on the combination K' + K of total momenta. Since all meson 
degrees of freedom are frozen, the operator has one-baryon and many-baryon 
pieces. Besides the standard nucleonic-current part there are additional parts 
involving the A-isobar. We take one-baryon and two-baryon contributions into 
account, shown in Figs. 2-4. The horizontal lines in the diagrams indicate 
that the meson exchanges are instantaneous. The exchanged mesons are the 
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Hn Ha 

Figure 1. Hilbert space considered. It consists of a purely nucleonic sector TCn and 
a sector Ha in which one nucleon is turned into a A-isobar, indicated by a thick 
vertical line. 



Figure 2. One- and two-baryon processes contained in the used e.m. current. The 
dashed horizontal line indicates instantaneous meson exchange. Diagonal ir- and p- 
exchanges are taken into account as well as the nondiagonal pivj and lott^j contribu- 
tions. In this figure only purely nucleonic processes are depicted. 



pion (jr) and the rho (p). The nondiagonal prrj and lot:"/ contributions, to de- 
noting the isoscalar vector meson omega, are also taken into account for the 
purely nucleonic current of Fig. 2. The current of Fig. 3 couples purely nucle- 
onic states with states containing one A-isobar; the hermitian adjoint pieces 
are not diagrammatically shown; the nondiagonal pirj and lottj contributions 
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to the two-baryon current involving one A-isobar are not considered. Among 
the contributions between A-isobar states only the one of one-baryon nature 
is kept as shown in Fig. 4. The current is derived by the extended S'-matrix 
method of refs. [12-15]; it satisfies current conservation with the correspond- 
ing 7r- and p-exchanges in the coupled-channel two-baryon interaction Hi of 
one-boson exchange nature. The current is systematically expanded up to first 
order in p/rriN, p being a characteristic baryon momentum and tojv the nucle- 
onic rest mass. 



Figure 3. One- and two-baryon processes contained in the used e.m. current. In this 
figure processes are depicted in which one nucleon is turned into a A-isobar. The 
hermitian adjoint processes are taken into account, but are not diagrammatically 
shown. 



Figure 4. One-baryon process contained in the used e.m. current. It is the only 
process between A-isobar states which is considered. 



In the perturbative spirit for the evolution of photo processes, the e.m. in- 
teraction Hf m - acts only once, whereas the hadronic interaction Hj has exactly 
to be taken into account up to all orders. We use hadronic channel states, seen 
in the initial and final states |iPj) and |/P/) of the photo reactions, in the 
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form 



|$ BPs ) = |B)|p B ), (2.5a) 
|$„(q)i/ a (JVd)K) = |0 a (q)i/ a (JVd))|K>, (2.5b) 



with the energies 



^s(ps) = m B c 2 + -^, (2.6a) 
bmjv 

q 2 K 2 

Sjvd(qK) = m N c 2 + m d c 2 + + ^— , (2.6b) 

m Ni m d an d being the rest masses of nucleon, dcuteron and trinucleon 
bound state. The internal trinucleon bound state is \B), which is normalized to 
1. The product nucleon-deuteron channel state in the three- nucleon cm. frame 
is \<p a (q)u a (Nd)} in the notation of ref. [1], v a denoting all discrete quantum 
numbers. In both cases the cm. motion is explicitly added to the internal 
motion; in the thrcc-nuclcon channel with a photon, the total momentum P is 
different from the total momentum K = p B of the three nucleons, bound in 
the trinucleon bound state \B). 

The matrix elements of the e.m. interaction require fully correlated hadronic 
states, i.e., 

|$ BPs ) =±iQG{E b (v b )±i{))\B)\ Vb ) : (2.7a) 

|*W (qK {Nd)K) = ±iOG(E Nd (qK) ± iO) (l + P) /V3|$ Q (qH (Nd)K) , 

(2.7b) 

with the full resolvent 

G(Z) - z- H ]-Hr < 2 ' 8 > 

the free Hamiltonian H a contains the motion of the center of mass; P sym- 
metrizes the nucleon-deuteron product state; the individual kinetic energy op- 
erators are of nonrelativistic form; they yield the eigenvalues of Eqs. (2.6). The 
hadronic state (2.7b) is normalized to the delta function without additional 
normalization factors. Since the hadronic interaction Hamiltonian Hi acts on 
relative coordinates only, the full resolvent reproduces the bound state \B) and 
correlates the nucleon-deuteron scattering state only in its internal part, i.e., 

\9&(q)v a (Nd)K) = ^(qK^d)) |K). (2.9) 
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2.2 S-Matrix for Trinuclcon Photo Reactions 

The kinematics of the considered processes is shown in Fig. 5. The figure also 
defines the employed notation for the individual particle momenta; PN,Pd,PB 




PN Pd PB k. 



Figure 5. Schematic description of the considered photo reactions of the deuteron. 
The lines for the two-baryon and three-baryon particles are drawn in a special form 
to indicate their compositeness. 



are derived from the zero components of those momenta, i.e., En(pn) = 
Pno c , Ed(pd) = Pdoc and Eb(pb) = Pboc; they are relativistic ones in contrast 
to those of the nonrelativistic model calculation in Eqs. (2.6) for which the 
same symbols are used; we shall make sure in the following that no confusion 
arises. 

The S'-matrix and the differential cross sections take the forms given in the 
following subsections. 

2.2.1 Radiative Capture in Nucleon-Deuteron Scattering 
We give various alternative forms for the S'-matrix elements: 

(fPflS^iPi) = (-i)(2^h)H(k 1+ p B -p N -p d ){s f \M rc \ Sl ) 

1 1 , 

x 77. ^g —^= 2.10a 

(27rfr)6 ^2fc 70 c 2E B (p B )2E N (p N )2E d {p d ) 

(fP f \S rc \iPi) = {-2m)8{k l0 c + E B (p B ) - E N (p N ) - E d (p d )) 

x (fP f \H^ m -iOG(E l +iO)\iP l ), (2.10b) 
(fP f \S rc \iP l ) = (-2ni)S(k 70 c + E B (p B ) - E N (p N ) - E d (p d )) 

x (B\ j"(-k 7 , p s + p N + p d )\ip { a +) (tn)v ai (Nd)} | qi _ p d -2 PAf ; 

(2.10c) 
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(sf\M rc \si) is the singularity-free matrix element for radiative capture, from 
which the differential cross section 



da rc = 



dLips(p N +p d ,k 7 ,p B ) 
Ac 2 sf(jp N ■ p d ) 2 - m 2 N m 2 c 4 



( s f\ M 7W7 g g T (2.11a) 



is obtained. Its dependence on the spin projections m Si and Mj i of nucleon and 
deuteron in the initial channel, collectively described by Si , and on the helicity 
A and on the spin projection M. Bj of photon and trinucleon bound state in the 
final channel, collectively described by s/, are explicitly indicated. (sf\M rc \si) 
is Lorentz-invariant and can be calculated in any frame. It is defined by 



(s f \M rc \ Si ) = ^(2nh) 3 ' 2 ^2E B (p B )2E N (p N )2E d (p d ) 

x (B|^(-k 7 ,p B +p A r + p d )e;(k 7 A)|^+)(qO^(A r rf)>L_P^P £ 

Hi— 3 

(2.11b) 

according to Eqs. (2.10). It contains our model description of the hadronic 
and e.m. interactions which will — for example with respect to experimental 
energies — not be fully precise. The model description uses nonrelativistic 
energies for all hadrons involved. Since the matrix element in (2.11b) can be 
calculated in any frame, we choose the cm. frame, i.e., pd = Pjv,Pb = k 7 , 
and the following computational strategy: 

• The experimental deuteron energy in the lab frame determines the total 
energy of the system and the relative momentum qj in the initial channel. 
This step is done using relativistic kinematics. 

• All hadron energies involved in the matrix element (2.11b) are chosen 
nonrelativistically and consistent with their model values. 

• Taking the computed trinucleon binding energy for calculating the rest 
mass m B and the average nucleon mass for the nucleon mass mjv, i.e., 
TOatC 2 = 938.919 McV, the magnitude of the photon momentum |k 7 | 
is determined employing the nonrelativistic forms (2.6) for the hadron 
energies, i.e., |k 7 c|+TOBC 2 +k 2 /6TOjv = TOjvc 2 +TOdC 2 +q|/(4TOjv/3). Since 
the model rest mass is not the experimental one, neither for 3 He nor 
for 3 H reactions, the photon momentum does not have the experimental 
value. 



The other building blocks for the differential cross section (2.11a) arc the 
Lorentz-invariant phase-space element dLips(pN + Pd 1 k 11 p B ) 1 i.e., 



dLips(p N + p d , k y ,p B ) 



d 3 k 7 d 3 p B 



(2^hf8{k 1 +p B -p N - Pd ) - 30 7 , y - - , (2.12) 

{2TTh) i 2k~ i oc{2iTh) i 2E B {p B ) 



and the incoming flux 4c 2 \J (p^ ■ pd) 2 — to^to^c 4 . In contrast to the matrix el- 
ement (sf\M rc \si) which contains our model assumptions, the kinematic parts 
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of the differential cross section, the Lorentz-invariant phase-space element and 
the incoming flux, have to be calculated with the true experimental and rela- 
tivistic energies of nucleon, deuteron and trinucleon bound state. In the cm. 
system the differential cross section takes the form 



da: 



d 2 k 7 



1 



1 



64tt 2 {hcf 



(sf\M rc \ Sl ) 



2 Ik. 



1 



IPJvl (E N (p N ) + E d (p N )Y 



(2.13a) 



That form of the cross section still describes the transition between specified 
spin states. The spin-averaged differential cross section is 



da: 



d 2 k 7 



M N .M dj \M Bf d2k 7 



(2.13b) 



2.2.2 Two-Body Photo Disintegration of the Trinucleon Bound State 

The definitions of the S-matrix (/P f\S \iPi), of the singularity-free matrix 
element (sf\M dls \si) and of the differential cross section are similar to those of 
radiative capture in Eqs. (2.10)-(2.13). However, the initial and final states are 
interchanged and e*(k 7 A) is replaced by e M (k 7 A). The singularity-free matrix 
element is defined by 



(s f \M dls \ Sl ) = ^(2irhf/ 2 ^/2E B (p B )2E N (p N )2E d (p d ) 

x (^ ) (q/)^ Q/ (iVd)|j^(k 7 ,p A r+p d + p B )e M (k 7 A)|B)| q/ = p i ^ m , (2.14a) 

from which the differential cross section 

do* = (s f \M^\ Si ) 2 ^ps(k. + PB PN , Pd ) ^ (2 ub) 

4C \PB ' fcy ) 

is obtained. Its dependence on the helicity A and on the spin projection M. Bi of 
photon and trinucleon bound state in the initial channel, collectively described 
by Si, and on the spin projections m Sf and Mj f of nucleon and deuteron in the 
final channel, collectively described by Sf, are explicitly indicated, (s f\M dt3 \si) 
is Lorentz-invariant and can be calculated in any frame. We choose the cm. 
frame and follow the computational strategy of Sect. 2.2.1 for radiative capture: 

• The experimental photon energy in the lab frame determines the total 
energy of the system and therefore the cm. momentum k 7 in the initial 
channel. This step is done using relativistic kinematics. 



• All hadron energies involved in the matrix element (2.14a) are chosen 
nonrelativistically and consistent with their model values. 
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• Taking the computed trinucleon binding energy for calculating the rest 
mass nig and the average nucleon mass for the nucleon mass mjv, i.e., 
mjq(? — 938.919 MeV, the relative momentum q/ of the final nucleon- 
deuteron system is determined employing the nonrelativistic forms (2.6) 
for the hadron energies, i.e., m^c 2 + to^c 2 + q 2 /(4mAr/3) = |k 7 c| + 
vhbc 2 + k 2 /6mjy. Since the model rest mass is not the experimental 
one, neither for 3 He nor for 3 H reactions, that nucleon-deuteron energy 
does not have the experimental value. 

The other building blocks for the differential cross section (2.14b) arc the 
Lorentz-invariant phase-space element dLips(k 1 + pb,Pn,P<i), i-e., 

dLips{k 1 + p B ,PN,Pd) = 

(27rfi) 4 <5(fc 7 +p B -p N -p d ) . ^fj^ , x 77. J^^t t — n! ( 2 - 15 ) 

(2nh) 3 2E N (p N ) (27rfi) 3 2^; d (p d ) 

and the incoming flux 4c 2 (ps-fc 7 ). In contrast to the matrix element {sf\M dls \si) 
which contains our model assumptions, the kinematic parts of the differential 
cross section, the Lorentz-invariant phase-space element and the incoming flux, 
have to be calculated with the true experimental and relativistic energies of nu- 
cleon, deuteron and trinucleon bound state. In the cm. system the differential 
cross section takes the form 



uu c.m. 



d 2 p 



N 



64tt 2 (He) 2 



'.Vl-V'"" «;> 2 ^T- — T2- (2.16a) 

l k 7l (fc 70 c + £ B (k 7 )) 



That form of the cross section still describes the transition between specified 
spin states. The spin-avcraged differential cross section is 



VN AM fl , M Nf M df HiV 



(2.16b) 



The S'-matrix of photo disintegration and radiative capture are related, 
due to time reversal, in corresponding kinematics situations |k 7 |c + E b (pb) = 
E N (p N ) + E d (p d ), i.e., \(fP f \S di '\iP i )\ = K/P/l^liPi)!, where the initial 
and final states of the matrix elements for photo disintegration and radiative 
capture are to be interchanged. That S'-matrix relation is derived in Appendix 
A. As a consequence, the differential cross section of photo disintegration and 
radiative capture are obtained from each other by 



d<J d Z. / da r c c m . _6\pn\ 2 ( 

d 2 p N I d 2 k 7 4 1 • '> 

2.2.3 Spin Observables 

Spin observables are defined as in ref. [2] . In this paper we consider polarization 
in the nucleon-deuteron system only. Thus, analyzing powers A in radiative 
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capture are 

Tr(M rc S a M rc A 

A=^- (2.18) 

TrlM rc M rc U 

The spin operators S a refer to the nucleon, i.e., S a — {<Ji} 7 or to the deuteron, 
i.e., S a = {Si, Sjk}] we take over the definitions of ref. [2]. For the deuteron 
analyzing powers the equivalent spherical tensor notation iTn and Tim will 
often be used as in ref. [2]. Details of the calculation of spin observables are 
also given in Appendix A. 
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3 Results 

The results are based on calculations derived from the coupled-channel two- 
baryon potential A2 of ref. [16]. The potential is an extension of the purely 
nucleonic Paris potential [17] extended in isospin-triplct partial waves for single 
A-isobar excitation. The A-isobar is considered to be a stable particle of spin 
and isospin 3/2 with a rest mass m^c 2 of 1232 MeV. The extension is almost 
phase-equivalent [18, 19] to the underlying nucleonic potential. The potential is 
used in the separably expanded form of ref. [1]. Partial waves up to total pair 
angular momentum 1 = 2 are taken into account. 

The AGS three-particle equations for the trinucleon bound state \B) arc 
solved as in ref. [1]. The resulting binding energy m^c 2 — 3itinc 2 is -7.702 MeV 
for the coupled-channel potential A2 and -7.373 MeV for the Paris potential, 
the purely nucleonic reference potential. If the Coulomb interaction is taken 
into account, as proper for 3 He, the binding energies shift to -7.033 MeV and 
-6.716 MeV; the e.m. shifts are calculated perturbatively in ref. [20], amounting 
to 669 keV for a coupled- channel potential with A-isobar excitation and to 657 
keV for the purely nucleonic reference potential. Nevertheless, we use the purely 
hadronic values -7.702 MeV and -7.373 MeV as model values for calculating the 
matrix elements (2.11b) and (2.14a) of the considered photo reactions, since 
the calculation leaves out e.m. effects in the hadronic wave functions. For 3 He, 
as needed in most described reactions, they differ slightly from the experimen- 
tal value -7.718 MeV. We shall describe proton-deuteron radiative capture for 
the two deuteron lab energies 19.8 MeV and 95 MeV; according to our calcu- 
lational strategy, the photon energies in the final states are 12.05 (11.72) MeV 
and 36.97 (36.65) MeV for the theoretical binding energy of 3 He with (without) 
A-isobar, the discrepancy with the values 12.07 MeV and 36.99 MeV arising 
from the experimental 3 He binding energy is quite small. The wave function of 
the bound state is expanded in terms of the (CS) coupled basis states defined 
in Appendix B.l; the expansion is done according to Appendix B.2; the infi- 
nite expansion is first truncated to include all states up to total pair angular 
momentum / = 6 in (Ij) coupling and then transformed to (CS) coupling. The 
convergence of the presented results with respect to that truncation is checked 
and found to be satisfactory for all reactions corresponding to deuteron lab en- 
ergies between 1 and 100 MeV: Additional wave function components usually 
change observables by less than 0.1% in the complete kinematic regime with 
the exception of sensitive spin observables with values close to zero: They may 
change within 1%. 

The AGS three-particle equations for the nucleon-deuteron scattering states 
\^a Z \q)iy a (Nd)} are solved as in ref. [4]; in contrast to ref. [2], ref. [4] uses real- 
axis integration for the solution. The wave functions are expanded in terms 
of the (Ij) coupled states defined in Appendix B.l; the terms of the infi- 
nite expansion according to Appendix B.3, contributing to the matrix element 
(ipa ■ ) (q/)^ Q , / (iV(i)|j Ai (k 7 ,pjv + p c i + p_B)e/j(k 7 A)|i?) of photo disintegration, arc 
limited with respect to total three-particle angular momentum J and with re- 
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spect to total pair angular momentum / by the angular momenta contained in 
the trinucleon bound-state wave function and in the considered current mul- 
tipoles. We note that also for radiative capture that matrix element of photo 
disintegration is calculated and time reversal is employed according to Ap- 
pendix A. 

The e.m. current is taken over from ref. [12]. Certain coupling constants, 
especially those referring to the A-isobar, are not up-to-date anymore; never- 
theless, we decided not to change them in the present paper in order to preserve 
a valuable historic continuity with our past calculations [12, 15, 21] for related 
observables. The technique for calculating multipole matrix elements is devel- 
oped in refs. [22,23]; a special stability problem arising in the calculation is 
discussed in Appendix B.5. Our calculation for photo processes with internal 
excitation energies of the three-nucleon system up to 50 MeV includes electric 
and magnetic dipole and quadrupole contributions, abbreviated in an obvious 
way by El, E2, Ml and M2. The magnetic multipoles are calculated according 
to Eq. (B.12b) of Appendix B.4 from the one- and two-baryon parts of the 
current. The electric multipoles are calculated in the Siegert form according to 
Eq. (B.13a) of Appendix B.4. The first term of Eq. (B.13a) is calculated from 
the one- and two-baryon parts of the current. The second term of Eq. (B.13a) 
is the physically much more important proper Siegert term; it is replaced by 
model energies and Coulomb multipoles; the contributing charge density has di- 
agonal single-nucleon and single-A contributions only; the nucleon-A transition 
contribution as well as two-baryon contributions are of relativistic order and 
are therefore omitted in the charge-density operator when calculating Coulomb 
multipoles. 

3.1 Calculational Tests 

3.1.1 Test of the Employed E.M. Current 

As test of the employed cm. current operator the magnetic form factors Fm 
of the trinucleon bound state are calculated. In fact, the obtained results are 
recalculations of our previous ones given in refs. [12, 15, 21] with updated tech- 
niques; they are shown in Figs. 6 and 7 individually for 3 He and 3 H and for 
their isospin-separated parts. Concerning the bound-state wave function, more 
mesh points and channels are included in the present calculation, significantly 
improving the results in and beyond the first diffraction minimum technically. 
Concerning the e.m. current, a sign mistake got detected in the two-baryon 
transition current from nucleonic to nucleon-A states; that contribution is in 
general very small, but becomes significant when all other contributions can- 
cel, i.e., in the region of the diffraction minimum. Both effects combined, the 
one arising from the improved bound-state wave function and the other arising 
from a corrected two-baryon transition current, yield changes in the magnetic 
form factors at momentum transfer Q > 4 fm -1 , decreasing the weight of the 
A-isobar for the magnetic form factors somehow, compared with the results of 
ref. [21]. 
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3.1.2 Test of the Calculational Apparatus for Radiative Capture 

A detailed comparison of the calculational procedure of this paper and of the 
one in refs. [7,8] is carried out. The comparison is done for a purely nucleonic 
calculation and with a heavily truncated potential; only the two partial waves 
1 Sq and 3 Si — 3 Di are kept and the separable expansion is simplified to include 
three and four terms only, i.e., it is of rank 3 and 4, respectively. With respect 
to the current, only the El multipole is employed in the long-wave length form 
Teiec mi (^7/^) according to Eq. (B.15b). Furthermore, the choice of momenta 
and energies for the current matrix element follows refs. [7, 8] in the rest of 
this subsection and not our own computational strategy as discussed in Sub- 
sect. 2.2. With the stated choice, the conceptional differences between the two 
calculations should be minimal and a detailed comparison is meaningful. On the 
other hand, the practical difference between the two calculations is enormous, 
especially with respect to the treatment of the permutation operator. 

The achieved agreement between the calculations of refs. [7, 8] and of this 
paper is quite satisfactory. We have compared the differential cross section and 
spin observables iTn, T20, T21, T22 and A yy . Fig. 8 shows characteristic results. 
For 10 MeV we show the case of poorest agreement, for other observables the 
results of both calculations are not distinguishable in a plot. For 45 MeV we also 
show the poorest case of agreement; the differential cross section and T22 have 
similar discrepancies, while the other observables are often not distinguishable 
in a plot. We conclude that our code for calculating trinucleon photo reactions 
is numerically quite reliable. 

3.1.3 Multipole Expansion of Current 

Our full calculations are based on a coupled-channel two-baryon potential with 
single A-isobar excitation and on a corresponding cm. current with the same 
coupling. The current is expanded into electric and magnetic multipolcs. The 
convergence of results due to that multipole expansion of the current is studied 
in this subsection. 

The first study is carried out for the full calculations of radiative capture 
with A-isobar excitation; the results for spin-averaged differential cross sections 
and for the spin observable A yy for 95 MeV, the highest considered deuteron 
lab energy, are shown in Fig. 9. The electric dipole operator of the Siegert form, 
i.e., the second term in Eq. (B.13a) for El, makes a significant contribution 
to most observables in the studied energy range. The first term of Eq. (B.13a) 
for El is unimportant for the differential cross section, but sizeable for spin 
observables; this observed fact is not documented in Fig. 9. The expansion with 
respect to the considered multipoles, i.e., multipoles up to E2 and M2, appears 
well converged. The contributions from the E3 and M3 multipoles, calculated 
tentatively, turn out to be very small for all observables, as is shown in Fig. 9. 

In Fig. 10 the multipole expansion is studied for the energy dependence of 
the two spin observables A yy and T20 at 90° lab scattering angle of the proton. 
The study is carried out for calculations with the purely nucleonic reference 
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potential. Up to 30 MeV deuteron lab energy the electric dipole operator El 
in the Siegert form, i.e., the second term of Eq. (B.13a), determines the ob- 
servables; in fact, its long- wavelength limit form Eq. (B.15b) suffices. Beyond 
30 MeV the first term in Eq. (B.13a) for El and the Ml and M2 multipoles 
become significant, but they tend to cancel each other; this fact, which we con- 
sider a non-understood coincidence, makes the result based on the El operator 
in Siegert form close to the full result. The contributions of E2 and higher 
electric or magnetic multipoles remain small. In comparison with ref. [8] we 
note in passing that beyond 30 MeV deuteron lab energy also the improved 
hadronic interaction with respect to high partial waves and additional ranks in 
the separable expansion gives rise to rather sizeable changes. 

3.1.4 Frame Independence of Matrix Elements (s f \M rc \si) and (s f \M dis \si) 
in Eqs. (2.11b) and (2.14a) 

The current matrix elements (s f\M rc \si) and (sf\M dts \si) are conceptionally 
Lorentz-invariant scalars and could therefore be calculated in any frame with 
identical results. We choose the cm. frame for calculations of this paper. How- 
ever, our calculational apparatus does not strictly respect Lorentz invariance, 
the results are frame-dependent. We therefore check that possible frame depen- 
dence of results and alternatively calculate the same current matrix elements 
also in the lab system. We do so for two-body photo disintegration at 36.9 MeV 
cm. photon energy, i.e., for the matrix elements (s f\M dls \si) . The recalcula- 
tion in the lab system uses the initial lab photon momentum k 7 and requires 
the redetermination of relative momentum of the final nucleon-deuteron 
system employing the nonrelativistic forms (2.6) for the hadron energies, i.e., 
mpfC 2 + iridc 2 + qj/(4mjv/3) + k 2 /6m,N = |k 7 c| + to^c 2 . The current operator 
needed for calculating the S'-matrix element (sf\M dls \si) is of nonrelativistic 
order; it has the same functional form as for the calculation in the cm. system, 
except for the single-baryon convection current. 

The El contribution to the matrix element (s f\M dls \si) is calculated ac- 
cording to Eq. (B.14a), the E2 contribution according to Eq. (B.12a) and Ml 
and M2 contributions according to Eq. (B.12b). Especially the two individ- 
ual terms of the El contribution show a frame dependence of the order 2%; 
however, that frame dependence cancels out for the sum of all multipoles con- 
sidered. Resulting observables of two-body photo disintegration at 36.9 MeV 
cm. energy are frame-independent up to a level of 0.01%. Thus, we conclude 
that all matrix elements (sf\M rc \si) and {sf\M dls \si}, calculated in this paper, 
are sufficiently frame-independent. 
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Figure 6. Magnetic form factors Fm of 3 He and 3 H as function of momentum transfer 
Q. The results for the interaction with A-isobar excitation are shown as solid lines, 
whereas the results for the purely nucleonic reference potential, the Paris potential, 
are shown as dashed lines. The experimental data are taken from ref. [24] 
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Figure 7. Isoscalar and isovector trinucleon magnetic form factors Fm as function 
of momentum transfer Q. The results for the interaction with A-isobar excitation are 
shown as solid lines, whereas the results for the purely nucleonic reference potential, 
the Paris potential, are shown as dashed lines. The experimental data are taken from 
ref. [24] 
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Figure 8. Comparison of the computational procedure of this paper with the one of 
refs. [7,8]. The comparison uses the two spin observables T22 and iTn of radiative 
capture at 10 MeV and 45 MeV deuteron lab energies, respectively. The scattering 
angle 8 cm is the cm. scattering angle of the photon with respect to the direction of the 
proton. The results use a simplified dynamics for the hadronic and e.m. interactions 
as described in Subsect. 3.1.2. The results based on the technique of this paper are 
shown as solid lines, the ones based on the technique of refs. [7, 8] as dashed lines. 
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Figure 9. Convergence of the multipole expansion of the e.m. current. The differential 
cross section and the spin observable A vv are sample observables. Contributions from 
multipole up to E3, M3 are presented for radiative capture at 95 MeV deuteron lab 
energy as function of the cm. scattering angle of the photon with respect to the 
direction of the deuteron; the El contributions include both terms of Eq. (B.13a). 
The results refer to full calculations with A-isobar excitation. 
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Figure 10. Spin observables of radiative capture at 90° lab scattering angle of the 
proton with respect to the direction of photon as function of deuteron lab energy. A yy 
is shown in the top, where the experimental data come from refs. [28] (open squares) 
and [27] (upper triangle); T20 is shown in the bottom, where the experimental data 
come from refs. [29] (open circles) and [30] (cross). The results refer to calculations 
with the purely nucleonic reference potential. For the solid curves all multipoles up to 
E2, M2 are included, for the dashed curves only El in the Siegert form of term two in 
Eq. (B.13a), for the dotted curves its long-wavelength limit according to Eq. (B.15b) 
and for the dashed-dotted curves El with both terms in Eq. (B.13a). In the top the 
solid and dotted curves are indistinguishable and are therefore not drawn separately. 
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3.2 Physics Results 

3.2.1 A- Isobar Effects in Radiative Capture 

This subsection presents full results for radiative capture. They are compared 
with those of a purely nucleonic reference calculation without A-isobar excita- 
tion which are derived from a purely nucleonic potential and a purely nucleonic 
current; the Paris potential is the nucleonic reference potential. 

First, in Figs. 11 and 12 we give sample results for 19.8 MeV and 95 MeV 
deuteron lab energies. In the figures we follow the standard convention for the 
presentation of data as introduced in refs. [25-27]: At 19.8 MeV the cm. photon 
angle refers to the direction of the proton, i.e., it is the cm. 3 He angle with 
respect to the deuteron direction; at 95 MeV the cm. photon angle refers to 
the direction of the deuteron. 

The total A-isobar effect is for all considered observables rather moderate. 
The inclusion of the A-isobar does not change the phase-space or flux contri- 
butions to the cross sections, they are purely nucleonic and are given by the 
experimental conditions. However, the inclusion of the A-isobar changes the 
current matrix element; its change can arise from four different sources: 

• The inclusion of the A-isobar changes the theoretical trinucleon binding 
energy and thereby changes the theoretical photon energy of the process. 

• The inclusion of the A-isobar adds A-components to the bound-state 
wave function and alters its nucleonic ones. 

• The inclusion of the A-isobar adds A-components to the wave function 
of the nucleon-deuteron scattering state and alters its nucleonic ones. 

• The inclusion of the A-isobar adds current components which couple to 
A-components in the hadronic wave functions. 

Figs. 13 and 14 try to distinguish the respective importance of those A-isobar ef- 
fects. The figures demonstrate that at the lower energy the A-isobar effects 
come mainly from changes in the binding energy and wave function of the 
bound state, at the higher energy mainly from the scattering state and the 
A-componcnts of the cm. current. 

Second, Fig. 15 shows results for selected observables, i.e., A yy and T20 at 
90° lab scattering angle of the proton, as function of deuteron lab energy, the 
A-isobar effect is small. We remind the reader that in Fig. 10 the multipole 
expansion is discussed for both observables. It is a physically remarkable effect, 
proven there, that below 30 MeV deuteron lab energy both observables are 
determined by the El Siegert operator of Eq. (B.15b) of the long-wavelength 
limit in a model-independent way. In fact, the contribution of that simple op- 
erator describes the whole energy regime considered extremely well, but that 
result is above 30 MeV deuteron lab energy due to a cancellation of rather 
sizeable other contributions arising from the first term of the El operator in 
Eq. (B.13a) and from other electric and magnetic multipoles. 
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In conclusion, we have theoretical predictions for existing experimental data 
in the whole considered energy domain. However, we refrain from showing more 
results; they are available from the authors upon request. 
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Figure 11. Radiative capture at 19.8 MeV deuteron lab energy as function of the 
cm. scattering angle of the photon with respect to the direction of the proton. The 
differential cross section is shown on the top; the experimental data are taken from 
ref. [25]. The spin observable T20 is shown on the bottom; the experimental data are 
taken from ref. [26]. The full results for the interaction with A- isobar excitation are 
shown as solid lines, whereas the results for the purely nucleonic reference potential, 
the Paris potential, are shown as dashed lines. 
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Figure 12. Radiative capture at 95 MeV deuteron lab energy as function of the 
cm. scattering angle of the photon with respect to the direction of the deuteron. The 
differential cross section is shown on the top. The spin observable A yy is shown on the 
bottom; both sets of experimental data are taken from ref. [27]. The full results for 
the interaction with A-isobar excitation are shown as solid lines, whereas the results 
for the purely nucleonic reference potential, the Paris potential, are shown as dashed 
lines. 




Figure 13. Contributions to the A-isobar effect for radiative capture at 19.8 MeV 
deuteron lab energy as function of the cm. scattering angle of the photon with respect 
to the direction of the proton. The cm. differential cross section is shown on the top; 
the spin observable T20 on the bottom. Shifts in the observables are shown. The full 
shift in the prediction arising from the different predictions of the nucleonic reference 
potential and the coupled-channel potential is shown as solid lines. That full shift is 
split up into four parts, i.e., the contribution arising from the theoretical change of the 
triton binding energy due to scaling (dotted lines), the contribution arising from the 
change of the nucleonic components of the bound wave function (dashed-dotted lines), 
the contribution arising from the change of the nucleonic components of the scattering 
function (short-dashed lines), and the contribution arising from the A-components 
in both baryonic wave functions and the e.m. current (long-dashed lines); those four 
contributions add up to the full shift making up the complete A-isobar effect. 
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Figure 14. Contributions to the A-isobar effect for radiative capture at 95 MeV 
deuteron lab energy as function of the cm. scattering angle of the photon with respect 
to the direction of the deuteron. The differential cross section is shown on the top; 
the spin observable A yy is shown on the bottom. Shifts in the observables are shown. 
The meaning of the different types of lines is the same as in Fig. 13. 
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Figure 15. Spin observables of radiative capture at 90° lab scattering angle of the 
proton with respect to the direction of photon as function of deuteron lab energy. A yy 
is shown in the top,where the experimental data come from refs. [28] (open squares) 
and [27] (upper triangle); T20 is shown in the bottom, where the experimental data 
come from refs. [29] (open circles) and [30] (cross). The results for the full interaction 
with A-isobar excitation are shown as solid lines, whereas the results for the purely 
nucleonic reference potential, the Paris potential, are shown as dashed lines. 
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3.2.2 Results for Two-Body Photo Disintegration 

According to Eq. (2.17) the cross sections for photo disintegration and radiative 
capture are simply related. In fact, the experimental data as shown in the figures 
of this subsection are taken from radiative capture; 12.2 (36.9) MeV photon 
cm. energy corresponds to 19.8 (95.0) MeV deuteron lab energy of radiative 
capture. 

The relation of Eq. (2.17) holds due to time reversal relation between ex- 
perimental cross sections. Of course, time reversal relation remains true also 
for the matrix elements M dls and M rc of the model calculations, however, the 
meaning of kinematic correspondence is different in model calculations. We 
voted to always determine the initial cm. momentum from the experimental 
lab energy by relativistic kinematics using experimental rest masses; the final 
cm. momentum and energy are then determined by energy conservation using 
the nonrclativistic kinematics of the model calculation. Thus, the kinematics of 
the final state does not match properly the experimental situation. Our recipe, 
which is created by the fact that the model calculation is based on nonrelativis- 
tic quantum mechanics and is unable to reproduce the 3 He binding precisely, 
clearly creates a conceptual uncertainty, the magnitude of which this subsection 
tries to estimate. 

We therefore calculated the spin- average differential cross section of 3 He 
two-body photo disintegration in two ways. In the first way, the experimen- 
tally corresponding radiative capture process is calculated and Eq. (2.17) is 
employed for the step to two-body photo disintegration. In the second way, 
two-body photo disintegration is calculated directly according to the proce- 
dure of Subsect. 2.2.2. The results are compared in Figs. 16 and 17 for a purely 
nucleonic reference calculation, since the purely nucleonic reference calculation 
misses the 3 He binding most and therefore should show the largest differences 
between the two calculational procedures. The complete A-isobar effect is also 
given in order to set the scale for the conceptual uncertainty of these calcu- 
lations. Results are shown for 12.2 MeV and 36.9 MeV photon cm. energies. 
Compared with the moderate complete A-isobar effect the conceptual uncer- 
tainty is small at the larger photon energy, but sizable at the small one. 

In Fig. 18 we show the energy dependence of the lab differential cross section 
for triton two-body photo disintegration with two-body break-up at 90° lab 
scattering angle of the nucleon with respect to the direction of the photon. 
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Figure 16. CM. differential cross section of 3 He photo disintegration with two-body 
break-up at 12.2 MeV photon cm. energy as function of the cm. scattering angle 
of the proton with respect to the direction of the photon, (a) In the upper part, the 
results of computational methods based on detailed balance according to Eq. (2.17) 
and on the procedure of Subsect. 2.2.2, are compared by the dashed line and the solid 
line, respectively. The results refer to a purely nucleon reference calculation, (b) In the 
lower part of the figure, the complete A-isobar effect is shown; here the calculations 
use the procedure of Subsect. 2.2.2. The full result for the coupled-channel interaction 
with A-isobar excitation is given by the solid line, the purely nucleonic reference result 
by the dashed line. 
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Figure 17. CM. differential cross section of 3 He photo disintegration with two-body 
break-up at 36.9 MeV photon cm. energy as function of the cm. scattering angle of 
the deuteron with respect to the direction of the photon, (a) In the upper part, the 
results of computational methods based on detailed balance according to Eq. (2.17) 
and on the procedure of Subsect. 2.2.2, are compared by the dashed line and the solid 
line, respectively. The results refer to a purely nucleon reference calculation, (b) In the 
lower part of the figure, the complete A-isobar effect is shown; here the calculations 
use the procedure of Subsect. 2.2.2. The full result for the coupled-channel interaction 
with A-isobar excitation is given by the solid line, the purely nucleonic reference result 
by the dashed line. 
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Figure 18. Lab differential cross section for 3 H photo disintegration with two-body 
break-up as function of the photon lab energy. The results refer to 90° nucleon lab an- 
gle with respect to the direction of the photon. The results of the full calculations with 
A-isobar excitation are shown as solid lines, those of the reference calculations with- 
out A-isobar excitation as dashed lines. Experimental data are taken from refs. [31] 
(open circles), [32] (open squares), and [33] (upper triangles). 
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4 Conclusion 

The paper calculates radiative nucleon-deuteron capture for selected energies 
and compares the theoretical predictions with existing proton-deuteron data. 
Some results for two-body photo disintegration are also given. The calculation 
includes the A-isobar in the hadronic and the e.m. interactions. It employs 
fully correlated hadronic states and an e.m. current consisting of one-baryon 
and two-baryon contributions. 

The A-isobar effect is isolated; it is moderate. At lower energies, it is mainly 
due to the shift in triton binding energy and due to the corresponding change 
in the nucleonic components of the bound-state wave function, brought about 
by the presence of the A-isobar. It is also found that below 30 MeV deuteron 
lab energy the spin observables A vv and T 2 o at 90° lab scattering angle can be 
obtained model-independently from the El Siegert term in the long-wavelength 
limit. At higher energies, the A-isobar effect arises dominantly from the changes 
in the scattering wave function and in the e.m. current due to the A-isobar. 
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Appendix A. Calculation of Amplitude 

In Subsects. 2.2.1 and 2.2.2 the calculation of cross sections for nuclcon-deuteron 
radiative capture and of two-body photo disintegration of the three-nucleon 
bound state is described. Both reactions are related by time reversal T. We 
use in Subsects. 2.2.1 and 2.2.2 the same notation for the momenta of both 
reactions, though they are in fact different, i.e., 



T\Vn) = 1 


(-Pn)) = 


lp* s ) , 


(A.la) 


t\p7) = 1 


(-P r /)) - 


\pf s ) , 


(A.lb) 


t\K) = 1 


(-k; c )> = 


|kf ) , 


(A.lc) 




(-PB )> = 


|pf s > • 


(A.ld) 



In this appendix we make the difference between the momenta of both reac- 
tions notationally explicit by the superscripts rc and dis. Time-reversal also 
changes angular momentum eigenstates \jm) by T\jm) = (— l) J+m |j(— m)}. 
Furthermore, the spin quantization axis, chosen in the experimental set-ups 
for the definition of spin observables, are different in both reactions. Fig. A.l 
illustrates the kinematic situations of both reactions related by time reversal. 

The current matrix clement of radiative capture is obtained from the corre- 
sponding one of disintegration; the transverse helicities of the real photon are 
easier implemented in the disintegration matrix element choosing the direction 
of the incoming photon as z-direction. 

(BM s |jn-k; c ,p^ + p^ + p^)e Al (k; C A)l^i +) (q rc )^M / m s (Ad)) 

x (^-\q L du )v a (-M I )(-m B )(Nd)\j*Qi* i ' ,pj} s +pf s +p*')e„(k^ X)\B(-M B )). 

(A.2) 

The derivation of Eq. (A.2) uses the general properties of current and 
photon- field operators under time reversal, i.e., 

rj^{x)r- 1 = J^Tx) , (A.3a) 

rA^{x)r- 1 = A^Tx) , (A.3b) 

and the condition 

^(-kA) = c„(-kA) . (A.3c) 

In contrast to Subsects. 2.2.1 and 2.2.2 the projections of total deuteron 
and trinucleon angular momenta M/ and Mb as well as the projection of the 
nucleon spin m s are notationally indicated in the hadronic states. All spin pro- 
jections in those states refer to the direction of the incoming photon momentum 
k* s in disintegration as quantization axis. 
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Figure A.l. Kinematic situations of nucleon-deuteron radiative capture and of two- 
body photo disintegration of the trinucleon bound state in the cm. system. The 
relation between momenta connected by time reversal and the respective beam axes 
are shown. 



The only spin observables studied in this paper are analyzing powers in 
radiative capture. In contrast to the quantization axis used for the matrix ele- 
ments (A. 2), the polarization of the deuteron and the nucleon are defined with 
respect to their respective momentum direction. Thus, either the spin depen- 
dence of the matrix element (sf\M rc \si) has to be turned to the new quantiza- 
tion axis or the deuteron and nucleon spin operators {Si, Sjk} and <7j have to be 
transformed to them by the rotation D{a(3^) = exp(—iaJ z /H)exp(—if3J y /h) 
x exp(—ijJ z /h), J being the general angular momentum operator creating the 
rotation, i.e., 

S l — > D{Qe cm Q)S l D\Qe cm Q) , (A.4a) 
Oi — » £>(0 - (tt - cm )O)<7i£>t (0 - (tt - cro )O) . (A.4b) 

The transformation of Sjk follows from Eq. (A. 4a). With those transformations 
the spin observables (2.18) of Subsect. 2.2.3 are calculated. When the depen- 
dence of observables on the photon helicity and on the trinucleon spin arc 
studied, no additional transformations of states or operators will be required. 
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Appendix B. Details of Calculation 

This appendix calculates the current matrix element of two-body photo dis- 
integration of the trinucleon bound state (ip~ (q/)^ Q/ (iVci) | j'' (k 7 , + Pd + 
p_B)e A1 (k 7 A)|B). The calculation is done in the cm. system, i.e., pd = — Pat, pb = 
— k 7 . The calculation uses a nonrelativistic potential model for the hadronic 
interaction. The baryon energies are employed in nonrelativistic forms. 

B.l Basis States 

The basis states \pqv) a describing the internal motion of three baryons are 
defined in Eq. (2.2) of ref. [1]. The discrete quantum numbers of the basis 
states, abbreviated by v are indicated in Fig. B.l. The subscript a at the states 
defines the role of particles as pair and spectator; this is why the subscript a 
at quantum numbers is omitted, unless ambiguities could arise. In Fig. B.l the 
(Ij) coupling scheme is used, i.e., 




Figure B.l. Three-baryon Jacobi momenta, and discrete quantum numbers. The 
spectator baryon is labeled a, the pair is made up of baryons (5 and 7. 

coupled with their respective spins S and s to total pair and spectator angular 
momenta / and j, which are then combined to total angular momentum J 
with projection A4j; the parity quantum number can be derived according to 
n = (— ) L (— The isospin coupling of the pair isospin T and of the spectator 
isospin t is done correspondingly for total isospin T with projection M.t- In the 
trinucleon bound states and in nucleon-deuteron scattering without Coulomb 
interaction only basis states with total isospin T = | are needed. The additional 
quantum numbers (Bb) give the baryon characteristics of the pair and of the 
spectator, B = 1(0) standing for a two-nucleon (nucleon-A) pair, b = \{—\) 
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for a spectator nucleon (A); baryon characteristics could be read off from the 
individual spin and isospin quantum numbers s and t; the additional quantum 
numbers (Bb) are introduced for convenience. The discrete quantum numbers, 
distinct from the continuous Jacobi momenta p and q, are abbreviated by 

Practical calculations suggest other coupling schemes for different parts of 
the numerics. The (£S) coupling scheme 

v{CS) = {(Ll)jC[(s s 1 )Ss}S}JMj{(t ( st 7 )Tt}TM T Bb (B.2) 
and the channel-spin coupling scheme 

v{lK) = {l([L(s l3 s 1 )S}Is)K}JM J {(t f ]t 7 )Tt}TM T Bb (B.3) 
are also used. The transformation between coupling schemes 

r l i c \ 

(v{CS)\v{Ij))=tSIjl S s S \, (B.4a) 
I I J J J 



{v{lKMIj)) = {-) l+K+I+i Kj{ I l j * }. (B-4b) 

(v{lK)(CS)\l/(CS)) = (-)L+S+s+K+S+J-LlgjK 

S L I I f L S K I (B.4c) 



K s S J I J I C 



using the abbreviation a = \J1a + 1, are unitary. They do not involve isospin 
and baryon content and arc independent from the continuous Jacobi momenta 
(jpq) at which the transformations are carried out. 

We calculate the matrix element (ip~(qf)i/ af (iV<i)|j A1 (k 7 , ps)e AJ (k 7 A)|i?) of 
photo disintegration first and then use time reversal for obtaining the corre- 
sponding one of radiative capture. The latter step is described in Appendix A. 
The matrix element has three ingredients, i.e., the current operator 
j M (k 7 , ps)e At (k 7 A), the initial bound state \B) and the final scattering state 
(ip~ (q/)f a/ (iV<i)|. Since the three baryons can be considered identical, whereas 
the employed basis \pqv) a is not totally antisymmetric, explicit antisymmctriza- 
tion has to be carried out with the permutation operator P = P123 + f*32i- The 
permutation operator P is described in ref. [2]. 

B.2 Three-Nucleon Bound State \B) 

The three-nucleon bound state \B) is obtained from the corresponding Faddeev 
amplitude \ip a ) according to 



\B)=N{l + P)\^ a ) 



(B.5a) 
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./V being the normalization factor ensuring (B\B) = 1. The Faddeev amplitude 
|V>a) is best calculated in (Ij) coupling, whereas for the computation of the 
current matrix element (£S) coupling is the better choice, i.e., 

( P qu(CS)\B) =N I ' dp'd q 'p /2 q' 2 J2{pq^jCS)\(l + P)\ P 'q'^(C'S')) 

J u'(C'S') I'j' 

x (v'(£'S')\v'(I'f)}(p'q'v'(rf)\ip a ). (B.5b) 



B.3 The Thrce-Nucleon Scattering State \ip^\q)v a {Nd)} 

The antisymmetrized nucleon-deuteron scattering state of internal motion 
|i/4 + ^(q)^a(-/VeQ) is obtained from the corresponding plane- wave state 
\4> a (q)v a (Nd)) by the full resolvent G(Z), i.e., 

|^i +) (q)^(iVd)) =iOG(E a (ci)+iO)(l + P)/V3\cl> a (ci)iy a (Nd)), (B.6a) 
|V>i +) (qH(^)) =(1 + -P)/V3[l + G (E + iO)T a (E + iO)G {E + iO)U(E + iO) 
x \M<l)v a (Nd))\ E=Ea{ci) . (B.6b) 

In Eq. (B.6b) the two-baryon transition matrix T a (Z) and three-particle multi- 
channel transition matrix U{Z) are introduced; they are defined in ref. [1]. In 
this section, as in ref. [1], the three-particle cm. part of the kinetic energy op- 
erator is left out; E a (q) is the internal nucleon-deuteron energy corresponding 
to the relative momentum q, i.e., E a (q) — ^^(qK) — K 2 /(6tojv) according 
to Eq. (2.6b). The two-baryon interaction is separably expanded. Thus, the 
plane- wave state has the following form 



= G (£ Q (q) + «0)| 3 i 4o " o7oTo )M 7 M To )|q So ™^o™ 4o 6o>a, (B.7a) 
|0 a (qK(Wd)) 

= G (E a (q)+iQ) Yl S \mx(Io.mjMjTM T ) a 

x (7oM / im j |JA4 i7 )(ToM To iom to |T MTHlmisom.^m^Y^ (q) (B.7b) 



with the quantum numbers 7To = +, Iq = 1, To = Mt = and so = to = &o = 
i. The choice of the spin projections Mj and m s specifies experimental condi- 
tions, and the choice of the isospin projection mt = +5 makes the calculation 
refer to proton-deuteron scattering. Eq.(B.7a) relates the deuteron state back to 

the corresponding form factor \ga° 7ToIoT °' > MjM To ) of the separable expansion; 
in Eq. (B.7b) the partial- wave coupled states |«o9x(7oj)nj7A^ jTA4r)ai intro- 
duced in Eq.(2.4) of Ref. [2] for the description of nucleon-deuteron scattering, 
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arc used. The practical computation of the scattering state 



\^ +) a (q)MNd)) = (1 + P)/V3G (E + iO) 



x [l + \ga)T a (E + iQ)(g a \G (E + iO)U(E + iO)G {E + iO) 



J E=E a (q) 



x E E \wx{hmJMjTM T ) 



TJJMjTMt jnijlmi 



x{IoMijmj\JMj){T M To t Q m to \T A4t) {Irrnsomsljm^Y^ (q) (B.8a) 
|^ +) (qK(JVd)> = (l + i , )/V3G (S + iO) E E 



EE / qldq c \l' c q c Xc{Ic]c)IiJMjTM T )(l' c \T{E + lQ,q c Xc{Ic]c))\lc) 



x (JoM/jmjlJA^j-) (T M To t m to \T M T ) (lm lSo m s \jra 3 )Y{ mi (q) (B.8b) 

proceeds according to the strategy suggested by the separable expansion: 

• The operator Go(E + iO)U (E + iO)Go(E + iO) is computed between states 
\iqx{Ij)^JMjTMr)a in (Ij) coupling; it and T a (Z) are diagonal in 
the three-particle quantum numbers (UJM jTMr) and independent of 
M j and M T ; for G (E + iO)U(E + iO)G a (E + iO) half-off-shell elements 
are required. Eq. (B.8b) uses the matrix elements 



{i c qcXc(Icic)\X njT (E + iQ)\iq X (I])) and (i' c \T (E + iO, q c Xc(Icjc))\ic) of 
G (E+iO)U(E + iO)G (E + iO) and of T a (E + iO) as introduced in ref. [2]. 



• The number of form factors |g Q ), i.e., the assumption on the dynamics, 
truncates the number of arising partial waves \iqx(I j)TiJ M. jTMr)a 
for each total angular momentum J. The truncation in J comes from 
the structure of the current, the matrix element of which is computed. 

• (US) coupled components a (pqv(£S)\ip^ (q)v a (Nd)) of the scattering 
state are determined; those components are most convenient for the full 
calculation of the current matrix element. The permutation operator P 
is applied in (Ij) coupling and the result is then transformed to (£S) 
coupling. 

The three-nucleon scattering state \ip^ (q)u a (Nd)) can be computed ac- 
cording to the steps explained in this subsection and the current matrix element 



TlJMjT Mr jmjlmi 



x \i qx(hj)njMjT M T ) a + 




DC 



x (icq c Xc(Ic]c)\X njT (E + iO)\iqx(Ij)) 



E=E a (q) 
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(■i/>i - ^(q/)^ Q/ (iV(i)|j Ai (k 7 , pjv+Pd+Ps) e ^(k 7 A)|i3) can be obtained from it and 
the bound state \B). We use that calculational avenue for numerical checks, but 
prefer a more direct alternative strategy for determing the current matrix ele- 
ment: 

( V4~ } (q/ ) v a , (Nd) | j " (k 7 , Pa t + p d + p B ) e M (k 7 A) | B) 

= (<t> a (q,)v a/ ME + iO)(l + P)/V3^(k 7 , p N + p d + p B )e A1 (k 7 A)|B) 



E=E af (q f ) 

(B.9a) 



fi(Z) = 1 + (7(Z)Go(Z)T Q (Z)G (Z), (B.9b) 

OO 

fi(Z) = ]T (PT a (Z)G Q (Z)) n . (B.9c) 



n=0 



Eqs. (B.9) are consistent with Eq. (B.6b) which relates the scattering state 
( c ii) v a i (Nd)) to the multi-channel transition matrix {/ (2? Q (q) + i0) . Since 
t/(Z) satisfies an integral equation [1], that property carries over to the current 
matrix element (B.9a), i.e., 

Sl(E + »0)(1 + P)/V3.^(k 7 , pjv + p d + P i?Mk 7 A)|B) 
= (1 + P)/V3j^(k 7 , p N + p d + p B )e^(k J X)\B) 
+ PT Q (£ + iO)G (E + iO) [o(£ + iO)(l + P)/V3j^(k 7 , Pa t + p d + p B )e AJ (k 7 A)|B)] 

(B.10) 

with = E af ((\f). The integral equation (B.10) has the same kernel PT a (E + 
«0)Go(i?-M0) as the corresponding one for U (E+iO) and can therefore be solved 
with the same nummerical techniques, just its driving term (l+P)/v / 3j M (k 7 , pn+ 
Pd + Ps)e^(k 7 A)|_B) is different. When taking the scalar product of the corre- 
lated state resulting from the integral equation (B.10) with the channel state 
(4> a (qf)i> af (Nd)\, the required current matrix element follows. This alterna- 
tive calculational scheme, based on (B.10), is the one usually employed for the 
numerics of the paper; it solves the integral equation for the scattering state 
\^a'\q)h' a (Nd)} implicitly when forming the current matrix element. 



B.4 Current Operator 

We calculate the current matrix element of two-body disintegration of the trin- 
ucleon bound states for the current operator j M (k 7 ,pAr + Pd + PB)f M (k 7 A) in 
the cm. system, i.e., Pat = —p dl Pb = — k 7 , assuming that the momentum of 
the real photon defines the z direction, i.e., k 7 = e 3 and A = ±1. The current 
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operator has the decomposition 



j M (fc 7 e 3 , -fc 7 e 3 )e M (fc 7 e 3 A) 



r 0) 

elecX 



(k-r/h) + \T^ agnX (k,/h) 



(B.11) 



into electric and magnetic multipoles, i.e., T^P (k) and Tml gn m Ak) with 
fc = fc 7 /ft and fc 7 = |k 7 | for the rest of this appendix. In the cm. system 
the formal parametric dependence on the total momentum K' + K = — k 7 is 
dropped. The multipoles are defined by 



^elec rrij (^) 



ll-LT^){k) + ^-J—T^^(k)\, (B.12a) 

(B.12b) 



± magn mjX^ ) ~ 1 raj W 

with 

i-iy 



(B.12c) 



In Eq.(B.12c) (k) is the spherical tensor corresponding to the spherical har- 
monics Yi mi (k),j^(hk) the spherical tensor corresponding to the spatial part 
of the corresponding current j^(fik,K.' + K). The electric multipole operator 
is employed in the Siegert form 



~ elec m 



(k) 



'2.7 + 1 



j hk 



(B.13a) 



with the Coulomb multipole 



coul m 



t (k) = T^\k) = { -^J d 2 kY m (k)p(hk), (B.13b) 



p(hk) being the Fourier transform of the charge-density operator with nor- 
malization p(Q) = 1. The form (B.13a) follows from current conservation; the 
hadronic Hamiltonian H includes the cm. motion of the baryons. As is well 
known, the Siegert form of the electric multipole operator takes exchange two- 
baryon currents into account, even if the Coulomb multipole is derived from a 
charge density, composed of one-baryon operators only as done in the calcula- 
tion of this paper. 

For many low-energy photo processes the electric dipole makes the dominant 
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contribution to the current. 



<V4- } (q/H, (Nd)\T { Jl mj (k,/h)\B) 
= (1>-(qj)v at (Nd)\(-)V3T£i\k,/h) - V2cT^ ul mj {k,/h)\B) (B.14a) 
= (^") (qy) Va/ (JVd) | (-)^TP; 11 (fey/n) 



That electric dipole operator is still exact. However, in applications it is often 
used in configuration space as in Eq.(B.14b) and in the long-wavelength limit, 

[2 11 

k-yx/h <C 1, in which the first term Tm] (kj/h) can be neglect and the regu- 
lar spherical Bescl function of order 1, ji(k 7 x/H) is approximated by fc 7 x/3S. 
Taking the charge operator p(x) as single-baryon operator 



the sum on i extending over all baryons, 73 (i) being the isospin projection with 
values ±1/2 for nucleons and with values ±3/2, ±1/2 for a A-isobar, e p being 
the positive elementary charge, the long-wavelength electric dipole operator 
takes the form 



{^-{c i} W s {Nd)\T^ m ) 1 {k 1 /K)\B) 

« -^^^-(q/K^Ar^l^l^ + TsW)^^). (B.15b) 



Our calculation is in momentum space. Thus, the use of the configuration space 
operator (B.15b) is unnecessary. Nevertheless, we are able to simulate the use 
of the approximate operator (B.15b) by a trick: 

• Calculate the matrix element (B.15b) for a fictions unrealistically small 
photon momentum /C 7 for which the approximation ji(K. 7 x/h) w JC^x/Sh 
with high precision. 

• Since the dependence on the physical photon momentum fc 7 is linear in 
Eq.(B.15b) the proper matrix element for the operator (B.15b) is obtained 
by scaling with the factor fc 7 //C 7 , i.e., 



B.5 Instability Problem in Calculation of Current Matrix Elements 

The physics results of this paper depend on matrix elements of the e.m. cur- 
rent J P (Q) according to Eqs. (2.4), (2.11b) and (2.14a). The computation of 




(B.14b) 




(B.15a) 



(^-(ci f )^Nd)\T^ c ll(k 7 /h)\B) 
= ^ (C (<1/K, (Nd)\T^ c mi (K,j/h)\B) . 



(B.15c) 
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those current matrix elements is quite involved and will not be spelt out in 
detail in this appendix. The computation is based on the multipole decomposi- 
tion of the current and on partial-wave decompositions of the initial and final 
baryonic states. The computation encounters particular numerical problems: 
The contributions arising from higher current multipoles and from higher an- 
gular momenta in the hadronic states get computationally split up into series 
of rather large numbers with alternating signs, the full results being compar- 
atively small; the straight-forward evaluation of those series becomes instable. 
We are unable to cope with that instability properly and to avoid it for good. 
Nevertheless, we believe that the results, presented in this paper, are compu- 
tationally quite reliable. This appendix tries to give convincing reasons for our 
claim. 

The current has one- and two-body contributions. For computational pur- 
poses, each current contribution is tensorially recoupled as tensor product of 
orbital, spin and isospin operators. The matrix elements of T^' a ' J '(fc) with 
k = k 7 /h take the general form 



(p'q'|r«^)(fc)|pq) 

= (~iy (4ir)i RK(a,b,e, d 2 k f p (k 7 p-,p+,v p ) f q (k,q-,q+,v q ) 

K •* 

(. - ( e )\ 0) 

F (K) (k)® \o (c \c p+ ,c p _,c pi c q )®S {d) \ J , (B.16) 
/ mi 



0( c ) being the orbital and being the spin part. The matrix elements of 
Tm] a ^ (k) are computed in partial 
three-baryon bound state \B), i.e. 



T^°' 3 ' (k) are computed in partial-wave basis with the operators acting on the 



(J'Mj,\T^\k)\B) 



1 



{-f'- M 'i ( J ' 3 2 ) {J'\\T^\k)\\B), (B.17) 



-Mj, rrij M 



B 



1/2 being the total angular momentum of the bound state and Mb its projec- 
tion. For a two-body contribution the formal structure of the resulting matrix 
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elements is as follows: 
(p'q'[{L'l')jC'(S's')S '}J'\ |T« J °W (jfe) | \B) 
= H)M47r)i X! ((^V)5'||^||(^)5) 

u={LlCSsS} 



xJ2 RK(a,b,e,j,I; K ) ^ RS(c,d,e,j;C',S',J',C,S,^;K,a) 

K <7 

^ ' RK{cp, Cp_ , Cp + , </9p, p) 




| p 2 ^(p'L'||p(>)( Cp+ ,p,^)|bL) 



oo 

X [ q 2 dq(q'l'\\Q^\c qi <p qi <p q )\\ql}(pqv\B). (B.18) 



X 





Both formulae (B.16) and (B.18) are admittedly intransparent ones. The quan- 
tities RK and RS are of geometric nature and just combinations of coefficients 
like 3j and 6j symbols arising from several recouplings. However, the detailed 
forms of (B.16) and (B.18) are irrelevant for the following numerical consider- 
ation. 

In Eq. (B.16) the complete momentum dependence is buried in the two 
functions f p (k,p-,p+,v p ) and f q (k,q-,q+,v q ). The function f p (k,p-,p+,v p ) 
is regular in k, the limit k = can be carried out without any problem, making 
that function independent of k and v p . In contrast, the function f q {k, q-,q+,v q ) 
carries a 5- function in k, which gets integrated out by acting the operator of Eq. 
(B.16) on the three-baryon bound state, i.e., by forming (p'q'|Tm^ a ' : ''(fc)|i?); 
then the limit k = can be performed leaving Y^(k) in (B.16) and in 
(p'q'|T^ ab) (fc)|B) the only dependence on k. As a consequence the tensor 
order 



k = 0, iffc = 0, (B.19) 

results. In Eq. (B.18) the characteristics of the considered current contribution 
enters by the non-negative integer parameters c p , c q and c with c q = or c q = 1. 
One numerically dangerous part of the expression is the ^-integration over the 
orbital operator Q^- lq \c q , tp q , tp q ) which weighs the bound-state wave function 
{pqv\B). Close inspection, e.g., according to Eq. (D.41) of ref. [23], shows that 
the chosen computation carries powers of the ratios q/k and q'/k which for 
momentum transfer k small compared with physically relevant q and q' in the 
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baryonic wave functions become unwantedly large. The numerical instability 
arises from the integration of those orbital operators Q^ Jq \c q , <p q , <p q ) for which 
the contributing tensor order j q becomes large. That source for instability is 
p-dependent and therefore contributes to the p-integration; it enhances insta- 
bilities in the integration over P^^ (c p+ , p,ip p ), which arise independently. We 
note, however, from the internal structure of the p-integral on P^ p \c p+ ,p,f p ) 
that in the limit k = the only contribution to P^ p \c p+ , p,f p ) is for vanishing 
ip p , i.e., 

ip p = 0, if k = . (B.20) 

The result follows from Eq. (D.14) of ref. [23]. Thus according to Eqs. (B.19) 
and (B.20) and according to the 9j-symbol in Eq. (B.18) ip q — as well as 
-f q = c q in the limit k = 0. Of course, outside the limit k = the tensor 
orders n > 0, tp p > 0, ip q > and j q ^= c q make contributions. Nevertheless, we 
observe a rapidly decreasing weight with increasing order 7 9 , i.e., we observe 
numerical convergence and stability well before instable fluctuations set in with 
high tensor orders 7 9 of Q^"\c q , <p q ,ip q ). 

In the table we give one example for the procedure and the reached nu- 
merical stability. The sample observable is the isovector magnetic form factor 
Fm(Q) of the three- nucleon bound state, the contribution arising from the ir- 
contact current is given. The observable in the table corresponds to the one 
plotted in the lower part of Fig. 7. The momentum transfer is called Q in that 
plot; it corresponds to the k in this appendix. The convergence with respect to 
the included maximal tensor order ~f q , m ax of Eq. (B.18) is studied. Numerical 
instability sets in when the tensor order j q . onS et is computed. 

Table: Contribution arising from the 7r-contact exchange current to the 
isovector trinucleon magnetic formfactor Fm{Q) 



Q [frn" 1 ] 





1 


, yq,max 

2 


3 


4 


'fq, onset 


0.001 


-0.622761 


-0.678250 


-0.678250 


-0.678250 


-0.678250 


6 


0.01 


-0.622743 


-0.678231 


-0.678231 


-0.678231 


-0.678231 


8 


0.1 


-0.620939 


-0.676272 


-0.676289 


-0.676289 


-0.676289 


> 8 


1.0 


-0.468388 


-0.510500 


-0.512039 


-0.511955 


-0.511955 


> 9 



Fig. B.2 gives another example for the onset of instability in the calcula- 
tion of typical multipole matrix elements needed for the description of photo 
processes. In the examples of the figure very high multipoles are chosen, since 
for them the instability problem becomes most severe. Though the examples 
are for special multipoles, special quantum numbers and for a final scattering 
state without final-state interaction, the on-set of instability is observed in the 
same general regime of 7 g for all matrix elements and for the required realistic 
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case in which final-state interaction in included. For lower multipoles the onset 
of instability occurs at somewhat higher j q . A corresponding on-set of insta- 
bility is seen, when calculating current matrix elements for inelastic electron 
scattering. 
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Figure B.2. Typical electric E3 (squares) and magnetic M2 (crosses) reduced mul- 
tipole elements in units of fm 3 ^ 2 , as function of j q ,max- The multipoles are derived 
from the two-nucleon 7r-contact current as in the table, i.e., they are isovectors with 
respect to their isospin character and they refer to radiative capture at 19.8 MeV 
deuteron lab energy. In the examples the matrix elements are calculated in the 
Wigner-Eckart reduced form with respect to total angular momentum and isospin 
for symmetrized plane-wave impulse approximation in the channel-spin coupling of 
Eq. (B.3), keeping the first term of Eq. (B.8b) only, i.e., without final-state inter- 
action: (i q X (lfK f )njT\\Go(E f + «))(1 + P)/ViT { Jl {magn) {k^h)\\B); the three- 
particle quantum numbers (JIJT) are chosen to be ((— ) 5/2 1/2), the spectator 
quantum numbers (IfKf) — (3 1/2); for easier graphical presentation the multipole 
matrix elements are scaled by the numerical factor 1.5 * 10 6 . 
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